Electric field driven donor-based charge qubits in semiconductors 
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We investigate theoretically donor-based charge qubit operation driven by external electric fields. 
The basic physics of the problem is presented by considering a single electron bound to a shallow- 
donor pair in GaAs: This system is closely related to the homopolar molecular ion Hj. In the case 
of Si, heteropolar configurations such as PSb + pairs are also considered. For both homopolar and 
heteropolar pairs, the multivalley conduction band structure of Si leads to short-period oscillations 
of the tunnel-coupling strength as a function of the inter-donor relative position. However, for any 
fixed donor configuration, the response of the bound electron to a uniform electric field in Si is 
qualitatively very similar to the GaAs case, with no valley quantum interference-related effects, 
leading to the conclusion that electric field driven coherent manipulation of donor-based charge 
qubits is feasible in semiconductors. 



Scalability constitutes one of the most attractive as- 
pects of solid state quantum computer (QC) propos- 
als. Group-IV or III-V semiconductor nanostructure- 
based QC architectures further benefit from the existence 
of a vast semiconductor microelectronics infrastructure. 
Semiconductor QC architectures typically rely on the 
electronic bound states of artificial atoms or molecules 
embedded in semiconductors, where the electron con- 
finement potential is provided by quantum dots (self- 
assembled or gate-defined) or shallow donor impurities. 
Information encoded in the electron spin or orbital state 
is processed through electronic charge manipulation by 
externally applied electric fields ^^^^2^ Single elec- 
tron coherent effects in GaAs have been demonstrated 
for electrons bound to donor impurities^ as well as elec- 
trons in double quantum dots^S*ii More recently, coher- 
ent time evolution of electrons in an isolated double quan- 
tum dot in Si was also reported^ 

Donor-based charge qubits are defined in close analogy 
with double dot charge qubitsjiS consisting of a singly 
ionized donor pair AB + & The electron confining poten- 
tial, usually taken to be parabolic in simple quantum 
dot models, becomes hydrogenic for donors. We have 
analyzed the tunnel coupling of P^" charge qubits in Si in 
Ref. \13L where it was shown that the multivalley struc- 
ture of the Si conduction band leads to an oscillatory 
and anisotropic behavior of the tunnel- coupling strength 
as a function of the donor pair relative position, favoring 
near-gap-closing cases. Nonetheless, charge qubit manip- 
ulation by electric fields performed in a Si double dotpi£ 
shows no indication of band-interference related phenom- 
ena in the tunnel coupling when compared to similar ex- 
periments performed in GaAs double dots^S It is thus 
important to clarify whether the experimental observa- 
tion indicates that in quantum dots the valley interfer- 
ence effect is absent, or that the electronic response to 
the external electric field is not sensitive to the valley 
interference, like the hyperfine coupling strength in Si:P 
system^4 



In this Communication we investigate the electric field 
manipulation of donor-based charge qubits in Si. For a 
better physical insight, it is instructive to consider ini- 
tially donor-pair molecular ions in GaAs. Although shal- 
low donor pairs in GaAs have simple and easy-to-control 
properties, sample preparation in GaAs poses an obvious 
difficulty that substitutional impurities at Ga or As sites 
will have entirely different behaviors (e.g. group-IV ele- 
ments act as donors when replacing Ga, but acceptors 
when replacing As). Therefore the GaAs results pre- 
sented here should be taken as a benchmark over which 
Si band effects are more clearly assessed. Furthermore, 
shallow donor binding energies in Si present a much wider 
distribution of values compared to GaAs, ranging from 43 
meV for Sb to 71 meV for Bi. Thus the optimal coupling 
of different donor species in Si (forming an heteropolar 
AB + molecular ion, which is feasible experimentally 15 ) 
would require the assistance of an external electric field, 
which again poses the interesting question on whether 
valley interference in these double donor systems might 
lead to any difficulty in the charge manipulation by the 
external field. Our goal here is to assess the prospect 
for the coherent manipulation of donor-based Si charge 
qubits, similar to that achieved in GaAsi2iii and Sii£ 
quantum dot based systems, in light of the knowni£ prob- 
lem associated with the quantum interference among the 
valleys in the Si donor-based systems. 

We model shallow donors in semiconductors within the 
effective mass approximation (EMA)^ where the elec- 
tronic properties of the host material are described by a 
few parameters, usually taken to be the static dielectric 
constant e and the effective masses at the conduction- 
band edge. The case of GaAs is particularly simple theo- 
retically, since the conduction band minimum at k = is 
non-degenerate and isotropic, characterized by an effec- 
tive mass m* = 0.067mo, where tuq is the free electron 
mass, and e — 12.56. We consider group-IV (VI) impuri- 
ties replacing the group-Ill (V) element Ga (As) , leading 
to a donor state, i.e., an additional electron in the system 
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FIG. 1: (Color online) Electric-field driven level anticross- 
ing for a single electron bound to a shallow donor pair in 
GaAs. The inset shows the calculated eigen-energies (the solid 
and dashed lines give the ground and first excited state re- 
spectively) versus applied axial field for interdonor distance 
R — 4.5 a*. As the field increases, the ground state changes 
from being spatially localized around donor A to B. Level an- 
ticrossing occurs at £ c = 0: In this situation the eigenstates 
are the symmetric and anti-symmetric superpositions |+) and 
|— ) (see text). The main graph gives the gap at anticrossing, 
A m i n , versus R. 

bound to a hydrogenic potential J7j(r) = — e 2 /(e|r — Rj|) 
centered at the impurity site R^. The one-electron 
Hamiltonian is written as H = i/caAs + Ui(r), where 
the first term describes the bulk semiconductor material 
(kinetic energy and periodic potential). Within EM A 
the electron bound to the donor is assumed to be highly 
delocalized in real space (with respect to the lattice pa- 
rameter), thus strongly localized in reciprocal space. Ad- 
ditional approximations 16 allow the electron wave func- 
tion to be written as ipi = Fj(r)</>k(r), the product of 
a slowly varying envelope function F by the band-edge 
Bloch state <^k which, in the case of GaAs, reduces to 
the rapidly varying periodic part of the Bloch function 
uq. The envelope function is the solution of the EM A 
Hamiltonian Hi — —h 2 y 2 /(2m*) + Ui(r), leading to 
Fi(r) = [l/V7r(a*) 3 ]exp(-|r - R;|/a*) with effective 
Bohr radius a* = e(ma/m*)ao « 10 nm and the eigen- 
energy for the donor state with respect to the bottom 
of the conduction band Ed = — [rn* / (e 2 mo)]Ky w —6 
meV, in excellent agreement with shallow donor binding 
energies m GaAs (see Tabled). 

The Hamiltonian for a single electron bound to an 
A — B donor pair at sites R,4 and H B in GaAs is 

H A _ B {£) = H GaAs + U A + U B + \e\£ s-r, (1) 

where the last term accounts for an external axial field 
E = £ s, with s = R/i? and R = — R^. For large 



TABLE I: Shallow donor parameters for GaAs and Si: Ed is 
the ground state eigen-energy with respect to the conduction 
band edge, and a* is the effective Bohr radius. 
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Sb 


P 


As 


Bi 


E D (meV) 


-5.8 c 


-43 


-45 


-54 


-71 


a* (nm) 


9.9 


1.32 


1.22 


0.97 


0.74 



"Effective mass approximation. 

''From Ref. Il7l where Ed is fitted to experiment and a* is calcu- 
lated within the multivalley spherical band approximation. 

c Experimental values of —Ed for C, Si, Ge, S and Se dono rs i n 
GaAs are, respectively, 5.9, 5.8, 5.9, 6.9 and 5.8 meV (from Ref. llSft . 



enough R (R > 4a*) this problem may be solved in anal- 
ogy with the LCAO solution for the Hj" molecular ion. 13 
The Hamiltonian is written in the non-orthogonal basis 
set {ipA, 4>b}-, which amounts to writing an EM A Hamil- 
tonian in the {Fa,Fb} basis set, since the periodic part 
of the Bloch function may be ignored in practice^^l 
All matrix elements are integrated analytically by sep- 
aration of variables using spheroidal coordinates*^ A = 
(ta + r B )/R and \i — (r A - r B )/R, with n = \r - Ri|, 
i = A, B. The electric field related term, not consid- 
ered in Ref. Il9l is transformed here through the relation 
§ • r = X^R/2. 

Figure ^ summarizes our results for donor charge 
qubits in GaAs. The inset gives the characteristic level 
anticrossing of the qubit states, driven here by the axial 
electric field, with the minimum gap A m ; n occuring at the 
anticrossing field £ c — 0. Far from the anticrossing re- 
gion, the eigenstates correspond to the electron localized 
around one of the individual donors, which we represent 
by \A) and \B) states, while for £ = £ c the ground state 
is the symmetric superposition |+) = (\A) + \B))/y/2 and 
the first excited state is the anti-symmetric superposition 
|— ) = (| A) — \B))/y/2. As the interdonor distance R in- 
creases, A min decreases smoothly, ranging typically from 
one to two orders of magnitude smaller than the single 
donor binding energy as R varies from 4a* to 8a*. 

The theoretical description of shallow donors in Si is 
more complex than in GaAs due to the six-fold degener- 
acy of the conduction band minimum of Si, corresponding 
to non-zero values of k^ t with fj, = 1, ... 6. Minima are lo- 
cated in reciprocal space 0.85 of the way along the T-X 
line. The energy dispersion near each of the band mini- 
mum is anisotropic, leading to different effective masses 
for the longitudinal and transverse directions along each 
of the T-X lines. The tunnel coupling of homopolar 
charge qubits in Si was investigated in detail in Ref. IT3I 
with the single donor wavefunctions taken in the Kohn- 
Luttinger (KL) EMA form, 16 i.e., for a donor at H A , 

Mr) = -7= E " RAK(r) e * k -( r - R -) , (2) 

where the envelope functions F^ir) are deformed IS* hy- 



FIG. 2: (Color online) Electric-field driven level anticrossing 
for the indicated A—B donor pairs forming an AB + molecular 
ion. The donors are positioned a distance R — 4.5 a P apart, 
along the [110] Si lattice direction, and the field is also applied 
along this direction. For the heteropolar cases, a positive field 
moves the electron away from the P site. 



drogenic orbitals. With respect to the GaAs results in 
Fig. ^ the most striking difference is the oscillatory and 
anisotropic behavior of A m i„ vs R for in Si. It is well 
establishedi2*22i2i that this behavior arises from interfer- 
ence between the plane- wave parts of the Bloch functions, 
which arc pinned to the respective donor sites. 

For the heteropolar donor pairs, the multivalley char- 
acter of the donor electron wavefunction is accounted for 
in the form of Eq. J2J). We adopt an additional simplifi- 
cation assuming isotropic envelope functions: F M (r) = 
[1/ y/ir(a* D ) 3 ] exp(-|r - K t \/a* D ), where D labels the 
donor species. We also assume that, as for the case of 
GaAs, the experimental donor binding energy satisfies 
Ed = {^d\Hd\^d), where He, is the single-donor mul- 
tivalley Hamiltonian. 13 Our treatment of the donor pair 
under an electric field follows closely the procedure out- 
lined above for GaAs: The donor pair Hamiltonian in 
the presence of an axial electric field, Ha-b(£), is writ- 
ten in the {i/ja^b} representation, in which one and 
two-center integrals are readily calculated in spheroidal 
coordinates. We adopt here the values of the effective 
Bohr radii calculated by Ning and Sahi£ for the 1S(A1) 
state of the different donors within the multivalley effec- 
tive mass spherical-band approximation. This approach 
involves a parametrization of the donor potential to fit 
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FIG. 3: (Color online) Calculated gap at anticrossing for 
biased AB + donor molecular ions in Si versus interdonor rel- 
ative position along the [100] (a), [110] (b) and [111] (c) lattice 
direction. Donor A is P and donor B is P (solid lines), Sb 
(dashed lines), As (dotted lines) or Bi (dot-dash lines). The 
axial electric field at anticrossing £ c , given in frame (d), is 
isotropic and depends on the donor species in the pair. Dis- 
tances are given in units of the P donor effective Bohr radius. 



the experimental values of Ed- In Table^we present the 
values of a* D and Ed adopted here. Comparison with the 
KL effective Bohr radii (2.5 nm and 1.4 nm) shows that 
a* D is probably underestimated. We expect our model to 
be more reliable regarding the energy estimates, since all 
donor binding energies are taken from experiment. 

Typical electric-field driven level anticrossing diagrams 
for donor molecular ions AB + are given in Fig. |2 where 
one of the donors (^4) is always chosen to be P. As com- 
pared to the Pj case, the eigen-energies for PSb + are 
higher while for PAs + and PBi + are lower, reflecting re- 
spectively the shallower or deeper nature of donor species 
B with respect to P. Since shallower donors have larger 
effective Bohr radii, yielding larger wavefunction overlap, 
the gap at anti-crossing A m i n also results to be larger for 
shallower donors, following the eigen-energies trend. We 
note that, as compared to the GaAs case given in the in- 
set of Fig. ^ no important qualitative difference emerges 
for the electric field dependence of the anti-crossing gap. 
Of course, for the heteropolar cases, level anticrossing 
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occurs at a non-zero bias field £ c . The energy scales for 
GaAs are typically one order of magnitude smaller than 
for Si, reflecting the difference in donor binding energies 
among these materials. 

Figure[3]gives the gap at anticrossing versus interdonor 
separation along three high-symmetry crystal directions 
in Si, as well as the corresponding bias field £ c . A ma- 
jor difference between these results and the GaAs case is 
the oscillatory and anisotropic behavior of A min vs R. 
Results for the molecular ion [solid lines in frames 
(a), (b) and (c)] are in excellent qualitative agreement 
with those reported in Ref. ^| (see Fig. 1 there) . ft is 
interesting to note that, even considering isotropic en- 
velope functions, anisotropic behavior is obtained due 
to interference among the plane-wave parts (pinned at 
sites Ha and Rg) of the six Bloch functions, corre- 
sponding to the wave vectors k M at the Si conduction- 
band minima. The results for heteropolar donor pairs, 
given by the dashed and dotted lines, follow the same 
overall behavior as the P^~ case, showing that the os- 
cillatory tunnel-coupling strength is inherent to the Si 
host, and persists in asymmetric molecular configura- 
tions. The same conclusions applicable to the Pj case^ 
in particular the tendency towards gap closing, remain 
valid here. Comparison of the solid line with the dashed 
and dotted lines in Fig.[3{a), (b) and (c) show that the 
coupling strength increases (decreases) when shallower 
(deeper) donors form the pair, as already noted in Fig. [21 
The bias field at anticrossing is plotted in Fig. Eld), 
where positive bias is arbitrarily chosen to lower the 
electronic potential at Re- Note that, for all heteropo- 
lar cases, \£ c \ decreases monotonically as R increases, 
which may be understood from the resonance condition 
(iPa\Ha-b(£c)\iPa) = (iPb\Ha-b(£c)\iPb)-- Each matrix 



element is taken between donor wavefunctions pinned to 
a single site, thus no quantum interference effects play a 
role, leading to the isotropic non-oscillatory behavior of 
£ C (R) in Fig.|2Id). The resonance condition in the lowest 
order approximation leads to \£ c \ — \Ea — Eb\/R, which 
fits the data presented in Fig. Etd) quite well. 

In this study we have shown that heteropolar donor 
pairs, like homopolar pairs, suffer from sensitive depen- 
dence of antocrossing gap on interdonor distance and ori- 
entation, resulting in a general narrowing of this gap, due 
to the valley interference in the Si conduction band. In 
the mean time, as summarized in Figs. and Eld), we 
have also shown that single electron response in a double 
donor system to an external electric field is not affected 
by the valley interference in the Si conduction band. In 
other words, applying an external field will not cause the 
tunnel coupling to exhibit oscillatory behavior as it does 
when inter-donor distance is varied. This conclusion is 
also applicable to Si double dot charge qubit, no matter 
whether its state structure is similar to that of a dou- 
ble donor in Si or double dot in GaAs. Given a fixed 
double-donor or double-dot geometry that has a large 
enough tunnel coupling, coherent charge manipulation 
in Si with electric field should be quite similar to that 
in GaAs, as recently demonstrated in Ref. ^3 via elec- 
trostatic pulse sequences in a Si double quantum dot. 
We conclude that external electric field induced coherent 
manipulation of donor-based charge qubits should be ex- 
perimentally possible in spite of the inter- valley quantum 
interference problemi&Si in silicon. 
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